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ABSTRACT 

We present observational evidence of apparent plasma rotational motions in the feet of a solar prominence. 
Our study is based on spectroscopic observations taken in the He I 1083.0 nm multiplet with the Tenerife 
Infrared Polarimeter attached to the German Vacuum Tower Telescope. We recorded a time sequence of spectra 
with 34 s cadence placing the slit of the spectrograph almost parallel to the solar limb and crossing two feet of 
an intermediate size, quiescent hedgerow prominence. The data show opposite Doppler shifts, ± 6 km s~\ at 
the edges of the prominence feet. We argue that these shifts may be interpreted as prominence plasma rotating 
counterclockwise around the vertical axis to the solar surface as viewed from above. The evolution of the 
prominence seen in EUV images taken with the Solar Dynamic Observatory provided us clues to interpret 
the results as swirling motions. Moreover, time-distance images taken far from the central wavelength show 
plasma structures moving parallel to the solar limb with velocities of about 10-15 km s~^ Finally, the shapes 
of the observed intensity profiles suggest the presence of, at least, two components at some locations at the 
edges of the prominence feet. One of them is typically Doppler shifted (up to ^ 20 km s~^) with respect to the 
other, thus suggesting the existence of supersonic counter- streaming flows along the line-of- sight. 

Subject headings: Sun: chromosphere — Sun: filaments, prominences 



1. INTRODUCTION 

Prominences are one of the most conspicuous solar struc- 
tures embedded in the ch romosphere or in the hot solar corona 
(for reviews see lLabrosse et al. 2010 and Ma ckay et al.ii20 10). 
They are made of cool magnetized plasma suspended above 
the solar surface. Recently, the understanding of promi- 
nence dynamics has greatly benefited from the high ca- 
dence and continuity of data taken with instruments such 
as the extreme ultraviolet light telescope (AIA; Lemen^eFal] 
20121) onboard NA SAs Solar Dynamics Observatory (SDO; 
Pesnell et al.l |2012|) . For instance, AIA observations show 
that the feet o f solar promine nces seem to be swirling in the 
lower corona (iLi et al.ir2012i) . Unfortunately, it is not pos- 
sible to infer plasma velocities from AIA narrow band im- 
ages. The only possibility is to measure local displacements 
of intensity structures in the plane of sky from consecutive 
images. This technique gives a measure of the plasma veloc- 
ity in the plane perpendicular to the line-of-sight (LOS), pro- 
vided the measured displacements can be considered as rea l 
plasma motions. It was recently used by Berger et alJ (1201 0|) 
to determine transverse velocities of dark bubbles in quies- 
cent filaments. Other authors have measured LOS plasma 
velocities from ground-based spectroscopic observations us- 
ingthe H^ line (e.g.,|Lrn et al. 2009, 2005; Zirker et al. 1998; 
Schmieder et al.1 119881: iMalherbe et al.. J983; Martres et al. 
198 Ih or lines in the EUV range (e.g., ^ngvold et al. 1985; 
Vial et al.lll979l:lLites et al][T976h . 

In this Letter, we measure spectral line Doppler shifts in the 
feet of an off-limb quiescent solar prominence using ground- 
based spectropolarimetric observations in the He 1 1083.0 nm 
triplet. We first present the observations and then determine 
Doppler shifts from time-averaged and individual profiles. 
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Fig. 1 . — SDO/AIA Fe XIV 21 1 A image of the quiescent hedgerow promi- 
nence. The straight line represents the TIP-II spectrograph slit, inclined 8.8° 
with respect the solar limb and at about 10^^ from it. The online journal 
mpeg animation shows the evolution of one of the prominence feet in the 
Fe IX 171 A band-pass during the TIP-II observations. 

We also determine plasma shifts parallel to the slit. We argue 
that these observations show evidence of swirling motions in 
the feet of the observed prominence. 

2. OBSERVATIONS 

The observations were taken between 12:53 UT and 13:58 
UT on 2011 May 19 with th e Tenerife Infrared Polarimeter 
(TIP-II;"Collados et al.'"2007') installed at the German Vacuum 
Tower Telescope (VTT) of the Observatorio del Teide (Tener- 
ife, Spain). TIP-II measured the four Stokes parameters of the 
He I 1083.0 nm multiplet in sit-and-stare mode. In particular, 
we placed the spectrograph slit (0'/5 wide and 80^' long with 
O.'^l? sampling) nearly parallel to the solar limb (forming an 
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Fig. 2. — Top panel: temporal variation (Y-axis) of the He I 1083.0 nm 
line peak intensity along the slit (X-axis). The feet of the prominence are 
the bright structures against a dark background. Bottom panel: time-slice 
diagram generated with the 171 A pass-band AIA data at the same location 
of the TIP-II slit. The feet of the prominence are seen as dark structures 
against a bright background in the AIA data. The resemblance between both 
maps suggests that AIA and TIP-II data were aligned successfully. 

angle of 8.8° with it). The exposure time per sHt position and 
polarization state was 7.5 s. The final cadence was 34 s and 
the total duration of the time sequence was 66.3 minutes (cor- 
responding to 1 17 slit images). The 1083.0 nm spectral region 
sampled by the spectrograph slit spanned 0.726 nm in 1.1 pm 
bins, thus including the Si 1 1082.7088 nm and an atmospheric 
water vapor line at 1083.21 nm. 

The TIP-II data reduction process includes dark current, 
flat-field, fringes correction, and polarimetric calibration. The 
absolute wavelength calibration was done using the atmo- 
spheric H2O line of the averaged spectrum as a reference. 
Then, we corrected the solar and terrestrial rotation and the 
relative shifts between the Sun and Earth orbits. This correc- 
tion varies 0.36 pm during the observations, introducing a sys- 
tematic error in the wavelength of ±100 m s~^ in 1 hr. Overall, 
the absolute error in the Doppler shifts is ~325 m s~^ 

In our observations, the TIP-II slit crossed two feet of a qui- 
escent hedgerow prominence located at the solar south east 
(Figure [T] sketches the position of the slit on a SDO/AIA 
Fe XIV 211 A snapshot). The position of the slit was kept 
rather stable during the observations thanks to the high perfor- 
mance of the Kiepenheuer-Instit ute Adaptive O ptics System 
installed in the VTT (von der Luhe et al.ll2003i) . The rms of 
the displacements of the slits in the plane of the sky should be 
less than 1 arcsecond. The evolution of one of the prominence 
feet can be seen in an Fe IX 171 A AIA movieQ, available in 
the electronic edition of this paper in the Astrophysical Jour- 
nal. The animation shows the apparent rotational motion of 
the prominence feet. 

^ In the movie, arrows represent the TIP-II slit. 



The temporal variation of the He I 1083.0 nm line peak in- 
tensity is displayed in the top panel of Figure |2l Note first 
that the He I 1083.0 nm triplet is seen in emission in off-limb 
prominences. The slit time sequence samples two feet of the 
prominence, the one on the right side is much more intense 
than the one on the left side. The data were obtained with rela- 
tively good seeing conditions. Indeed, the feet show a glimpse 
of the fine- scale in the prominence. We estimate a spatial res- 
olution of about I" using cross-correlation techniques. 

The bottom panel of Figure [2] displays a time-slice diagram 
generated with the SDO/AIA Fe IX 171 A data at the same 
location of the TIP-II slit. The position of the slit with respect 
to the AIA data was estimated using the orientation of the slit 
with respect the plane of the zero solar meridian, the angle 
between the slit and the line tangent to the nearest limb posi- 
tion, and Ha slit-jaw images. The co-alignment error should 
be within the AIA resolution limit of 1 .''5. In the AIA data the 
feet of the prominence are seen in absorption mainly from 
neutral hydrogen and neutral and ionized helium against a 
bright background. The He I 1083.0 nm peak intensity and 
the AIA absorption show similar spatial and temporal varia- 
tions. The reason for the good resemblance between the EUV 
and infrared data is that the EUV spectrum below 504 A ion- 
izes helium which recombines again producing an enhance- 
ment of the populations in its triplet states and in particular in 
the metastable level (2s^Si) of the transition. This produces 
an increase in the emission of the He I 1083.0 nm multiplet, 
which would be hardly explainable with collisional excitation 
alone (lAndretta & Joneslll997l:ICenteno et al.ll2008h . 

We used STEREO-B EUVI 195 A observations 
(iKaiser et al.l |^08) to estimate the filament orientation 
(prominences are seen in absorption when seen on disk). The 
prominence spine was rotated '^17° counterclockwise with 
respect to the solar south-north direction. Stereo data also 
helped to determine the solar longitude of the prominence, 
about 80° East. 

3. CALCULATION OF DOPPLER SHIFTS IN THE PROMINENCE 

FEET 

3.1. Average Doppler shifts 

We first restrict the analysis to the temporally averaged in- 
tensity profiles. Figure [3] shows the spatial variation of the 
peak intensity of the temporally averaged spectra along the 
slit. The intensity variation exhibits two large bumps corre- 
sponding to the prominence feet (one from 15'' to 30'' and 
another from 45" to 64"). Note that the peak intensity of the 
He 1 1083.0 nm triplet never reaches zero in these data. Hence, 
we constrained our analysis to those pixels whose peak inten- 
sity surpasses 30% of the maximum signal (dashed line in top 
panel from Figure O. The Doppler shifts are derived from the 
intensity profile by fitting a Gaussian. We fit three parame- 
ters: the Doppler shift, the amplitude of the Gaussian, and its 
width. 

The bottom panel of Figure [3] shows the inferred Doppler 
shifts, where positive means redshifts. If we focus on the se- 
lected pixels only, we readily notice that the left sides of the 
feet are blue-shifted whereas the right sides are red-shifted. 
The inferred Doppler shift smoothly varies from one edge to 
the other edge, being at rest at about the center of the feet. The 
red-shifts are greater in the right feet, around slit positions 
~60". There, the He I 1083.0 nm peak intensity decreases 
rapidly and the peak signal falls below the 30% threshold at 
64", coinciding with a dramatic decrease in the Doppler shift. 



First detection of apparent swirling motions in the feet of solar prominences 
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Fig. 3. — Top panel: variation of the time-averaged peak intensity of the 
He 1 1083.0 nm line along the slit, normalized to its maximum value. Bottom 
panel: Doppler shifts inferred from the He I triplet. In both panels, the vertical 
lines delimit the feet of the prominence (see text). 

3.2. Doppler analysis of individual profiles 

The visual inspection of the individual emission profiles 
corresponding to the prominence feet revealed that they show, 
at least, two different components. The one showing the 
greatest peak intensity (hereafter dominant component) is 
generally slightly blueshifted or redshifted depending on the 
location of the pixel with respect to the center of the promi- 
nence foot. The other (hereafter weaker) component show- 
ing a lesser peak intensity is conspicuously Doppler shifted. 
Examples of these profiles are displayed in Figure (H Their 
location is marked in Figures [2] and \5\ Profile (a) represents 
a prototypical, one-component He I 1083.0 nm emission pro- 
file, i.e., showing two peaks corresponding to the red and blue 
components of the He I fine- structure. Profile (b), located in 
the right side of the right foot, shows two components. In 
this case, the dominant component moves away from the ob- 
server in agreement with the Doppler pattern derived from the 
temporally averaged profiles while the other component is, 
surprisingly, strongly blue-shifted; i.e., moving towards the 
observer at about 12 km s~^. Not all profiles in that part of 
the prominence show the same behavior but there are few ex- 
ceptions in isolated places. For instance, the dominant com- 
ponent in profile (c) shows an opposite Doppler shift to that 
inferred from the temporally averaged profiles. Profile (d) 
also shows two components and strong Doppler shifts. In this 
case the profile belongs to the left side of the right foot; i.e., 
where the plasma moves toward the observer. Here, the dom- 
inant component is blue- shifted while the weaker component 
is red-shifted. In this example the Doppler shifts of the domi- 
nant component are quite relevant, up to ~ 10 km s~^ . Profile 
(e), belonging to the left foot, shows a similar behavior as pro- 
file (b) although here the shift of the dominant component is 
less conspicuous, probably because the intensities in the left 
foot are overall weaker than the right foot, while the weaker 
component shows a strong Doppler shift, of ^ -14 km s~^ 
Finally, the dominant component of profile (f), corresponding 
to the left side of the left feet, shows Doppler shifts opposite 
to that inferred from the temporally averaged profiles. 



The Doppler shifts have been determined by a two- 
components Gaussian fit to the observed spectral lines. The 
fits can be seen superposed on the profiles shown in Figured 
The fits are good although sometimes the blue component of 
the He I triplet is miss-fitted, probably because of the simplic- 
ity of the model approach. To model the fine- structure of the 
He I triplet we assumed that each of the two components con- 
sisted of three individual Gaussians with the same Doppler 
width, one for the blue component of the line and two for 
the red component of the line. For simplicity, we fix the ra- 
tio between the blue and the red components of the triplet to 
5. However, this ratio may change depending on the amount 
of ionizing EU V coronal irradiation and on the line opacity 
(ICenteno et al.i r2008), hence producing the misfit shown in 
the example profiles of Figure \5\ To select pixels with two 
plasma components we set a threshold in the full width at half 
maximum of the observed line profiles because profiles with 
two components tend to be much wider than the rest. This 
threshold is completely arbitrary but direct inspection of the 
data shows that most two-component profiles were selected. 

The inferred Doppler shifts from the two-component Gaus- 
sian fit are displayed in Figure \5\ The Doppler shift pattern 
of the dominant component (top-left panel) shows that the 
two sides of the feet move in opposite direction. This be- 
havior persists in time, being much stronger at the end of the 
observations (see top of upper left panel in Figure [5]) where 
the Doppler shift reaches ±6 km s~^ . The Doppler shift pat- 
tern matches the one obtained from the temporally averaged 
profiles. The Doppler shift map corresponding to the weaker 
component (top-right panel) shows red- shifts on the left side 
of the feet and blue- shifts when they are found on the right 
side of the feet, although exceptions exist, e.g., areas around 
pixels (c) and (f). The Doppler shifts may be as large as 
20 km s~^ . Interestingly, the Doppler shifts exhibited by these 
weaker components are of opposite sign to that obtained from 
the dominant component. 

In Figure [5] we also display the Doppler width of the domi- 
nant component inferred from the two-component Gaussian 
fit. The Doppler widths are about a factor of 1.5 larger 
than that expected from thermal broadening only, about 0.023 
nm at lOOOOK. Thus the observed emission profiles may be 
broadened because of the presence of unresolved components, 
either along the LOS or within the pixel resolution. Here we 
distinguish, at least, two components. Figure[5]also shows the 
optical thickness of the He I 1083.0 nm emission line. It has 
been estimated from a on e component fit to the data using the 
HAZEL inversion code (lAsensio Ramos et al.|[2008) . They 
are of the order of 0.2^0.6 around the position of the promi- 
nence feet. This means that the He I 1083.0 nm line is op- 
tically thin and therefore we cannot discard that the Doppler 
shifts we are measuring might be associated with real plasma 
velocities along the LOS. 

3.3. Transverse shifts 

Figure \5\ (bottom panels) represent time-distance maps 
taken in the blue and red side around the He I line cen- 
ter (1083.029 nm) and corresponding to Doppler shifts of 
±17 km s~^. To enhance the temporal variation of the struc- 
tures we have removed the averaged intensity profiles along 
the slit in the two maps. Both maps show intensity patterns 
moving along the slit. These patterns, often associated with 
flows, allow us to infer transverse shifts, i.e., plasma displace- 
ments parallel to the solar limb and perpendicular to the LOS. 
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Fig. 4. — Example He I 1083.0 nm individual profiles taken from different regions of the prominence. Panel (a) represents a prototypical example of a 
rest profile, showing the fine structure of the multiplet, with a weak blue component at 1082.909 nm (^Si--^Po) separated about 0.12 nm from the other two 
components at about 1083.029 nm (^Si-^Pi and ^Si--^P2), which are indistinguishable. Panels (b), (c), (d), (e), and (f) display individual profiles showing, at 
least, two-components. Dots represent the observations and solid a two-component Gaussian fit. Arrows in Figures [2] and |4] show their location with respect 
the prominence feet and the corresponding observing time. The dot-dashed and dashed vertical lines stand for the wavelength location of the shifted and main 
components. The dotted lines show the rest position. 



Two clear examples of such patterns are pinpointed by arrows 
with solid arrowheads in Figure [51 The slope of the structures 
implies horizontal speeds of ~10 km s~^ and 15 km s~^ . Inter- 
estingly, the profiles associated with these structures, profiles 
(e) and (f), are those where the dominant component showed 
a Doppler shift opposite to that obtained from the analysis of 
the temporally averaged spectra. 

4. DISCUSSION 

The Doppler shifts inferred from the temporally averaged 
intensity profiles show that the left and right sides of the 
prominence feet move in opposite directions, with a smooth 
variation from -2 km s~^ at the left edge to 4 km s~^ at the 
right edge in the right foot and from -1 km s~^ at the left 
edge to 1 km s~^ at the right edge in the left foot. We also 
inferred Doppler shifts from individual profiles using a two- 
components Gaussian fit. The results for the dominant com- 
ponent show speeds of about ±6 km s~^ at the edges of the 
prominence feet, again with a smooth variation along the slit. 
The sign of the Doppler shifts corresponding to the dominant 
component coincides with that obtained from the temporally 
averaged data. The inferred Doppler pattern is similar in both 
prominence feet and the pattern persisted during the whole 
observation in the right feet. 

To determine whether the measured Doppler shifts can 
be ascribed to real plasma velocities we inferred the opti- 
cal thickness of the observed spectral lines and their Doppler 
widths. The optical thickness reveals that the observed spec- 
tral lines are optically thin while the inferred Doppler widths 
suggest that there might be unresolved components within the 
spatial resolution limit and/or along the LOS. In these data 
two of these components are clearly distinguishable in some 
areas at the prominence feet. Hence, we believe that the in- 
ferred Doppler shifts are not contaminated by the superimpo- 
sition of fine structure elements. Therefore, we assume that 
the observed He I 1083.0 nm triplet is tracing real plasma 
velocities. The fact that the time evolution and spatial pat- 



tern showed by the EUV and He I 1083.0 nm data are similar 
strengthens this assumption. 

Since in the EUV movie one sees plasma apparently 
swirling around the prominence feet axis we argue that 
the inferred Doppler shifts are consistent with the presence 
of prominence plasma structures rotating counterclockwise 
around the vertical axis to the solar surface as seen from 
above. The inferred Doppler shifts along the slit would corre- 
spond to the LOS component of such rotational motion. The 
rotation is maintained during the whole observation. If the 
prominence feet are cylindrical, vertical structures with an av- 
erage radius of 20", it would take about 4 hr to complete a 
rotation at an angular speed of 6 km s~^ . This rotation speeds 
coincide with those i nferred in so-called solar tornados us- 
ing indirect methods (iLi et al.ll2012|; IWedemeyer-Bohm et al.l 
I2012HSU et ani2012l) . We also detect plasma shifts parallel to 
the slit using time- slice diagrams. 

Finally, the results show the presence of strong LOS flows 
(associated with the weaker component of the observed pro- 
files) localized in areas surrounding the prominence feet. The 
Doppler shifts reach up to ±20 km s~^ and have sign oppo- 
site to the Doppler shifts determined from the dominant com- 
ponent of the profiles. Interestingly, these flows are close 
to supersonic. The sound speed in prominence plasma tak- 
ing into account partial ionization and electron densities in 
the range from 10^ to 10^^ cm~^ and temperatures between 
6000 and 10000 K ran ges from 8.5 km s"^ to 14 km s"^ 
(c.f., iBerger et"an 120101) . The inferred Doppler shifts repre- 
sent lower limits since they correspond to the LOS component 
only. We have no clues about the origin of such strong LOS 
flows, although, the most probable s cenario is the existence of 
supersonic counter- streaming flows (|Zirker et al.lil99 8) along 
the line-of- sight in the prominence feet. The physical origin 
of the counter- streaming flows may be the Rayleigh-Taylor in- 
stability recently found in prominences (Berger et al. 20®. 
In this case, we would be detecting a combination of hot 
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Fig. 5. — Top panels: Doppler shifts resulting from the two-components Gaussian fit of the intensity profiles. Left panel represents the shifts of the dominant 
component, saturated at ±6 km s~^ whereas the right panel shows that corresponding to the weaker component, saturated at ±20 km s~^ (see Sect l3.2b . White 
areas in the top-right panel represent locations where the profiles show one dominant component only. Left middle panel displays the Doppler width of the 
dominant component retrieved from the two Gaussian component fit. Right middle panel represents the optical depth of the He I 1083.0 nm multiplet inferred 
from a one component inversion with HAZEL. Bottom panels: time-distance maps taken at zb 17 km s~^ (left and right, respectively for -i- and -) from the He I 
line center (1083.029 nm). Vertical lines delimit normalized intensity values larger than 0.3 units (see Figure[3). Arrows show the location of profiles displayed 
in Figurelll Arrows with solid arrowheads show locations of transverse Doppler shift patterns (Section 3.3). 



(upward) and cold (downward) plasma flows resulting from 
the Rayleigh-Taylor instability. Since these instabilities have 
been found acting in the direction perpendicular to the solar 
surface, projection effects may be also playing a role. 
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Ruiz Cobo and A. Diaz for constructive criticisms of the 
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